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O Abstract 

o 

Tj" ' In neuroscience, optics and condensed matter there is ample phys- 

. ical evidence for multistable dynamical systems, that is, systems with 

Q>^ I a large number of attractors. The known mathematical mechanisms 

that lead to multiple attractors are homoclinic tangencies and stabi- 
lization, by small perturbations or by coupling, of systems possessing 
I a large number of unstable invariant sets. A short review of the ex- 

Q ' istent results is presented, as well as two new results concerning the 

^ ■ existence of a large number of stable periodic orbits in a perturbed 

Q ! marginally stable dissipative map and an infinite number of such orbits 

in two coupled quadratic maps working on the Feigenbaum accumu- 
lation point. 



1 Multistability. Motivations. 

That are three basic motivations to study multistable systems, that is sys- 
tems that possess a large number of coexisting attractors for a fixed set of 
parameters. First there is ample evidence for such phenomena in the natural 
sciences, with examples coming, among others, from neurosciences and neu- 
ral dynamics @ § [|] § 1 § , optics @ g, chemistry HI p 
condensed matter [|T3| and geophysics |14 . 
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The second motivation to study multistable dynamical systems is the 
mathematical challenge of identifying the universal mechanisms that lead to 
multistability and to prove rigorously under what circumstances the phe- 
nomenon ,may occur. 

The third and last motivation arises from the field of control and tech- 



nological design. After the pioneering work of Ott, Grebogi and Yorke [15 



the field of control of chaos became a whole industry (|]T6l and references 
therein). Control of chaos deals with the local control of unstable orbits, 
either to achieve their stabilization or, alternatively, to target a dynamical 
system to some desired final state. Reliable stabilization of unstable peri- 
odic orbit requires either the knowledge of a good model of the system or an 
accurate local reconstruction of the dynamics. This is feasible in some low- 
dimensional systems, but it seems rather problematic for high-dimensional 
ones. The typical situation in control of chaos, is that of a strange attractor 
with an infinite number of embedded periodic orbits, all of them unstable. If, 
instead of an infinite number of unstable periodic orbits, one has, for exam- 
ple, an infinite number of sinks, the controlling situation would seem more 
promising. The sinks would of course have very small basins of attraction. 
Nevertheless, the control need not be so accurate, because it suffices to keep 
the system inside the desired basin of attraction. This, in principle, makes 
for a more robust control. 

In this paper I will concentrate mostly on the rigorous mathematical re- 
sults that concern multistable dynamical systems. Sects. 2 and 3 contain 
brief reviews of the Newhouse phenomenon and of the creation of a large 
number of periodic orbits by dissipative perturbations of conservative sys- 
tems. Finally, in Sects. 4 and 5, some new results are presented which, in 
Sect 3, use the techniques of deformation stability and, in Sect. 5, prove the 
existence of an infinite number of sinks for two coupled quadratic maps. 



2 DifFeomorphisms with homoclinic tangen- 
cies.The Newhouse phenomenon and be- 
yond 

Contrary to earlier conjectures that generic systems might have only finitely 
many attractors, Newhouse fl^ |0 proved that a class of diffeomor- 
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phisms in a two-dimensional manifold has infinitely many attracting peri- 



odic orbits (sinks), a result that was later extended to higher dimensions ||20 
Concretely, for two-dimensional manifolds the result is: 

Theorem (Newhouse, Robinson ||21||) Let be a map in a 2-dimensional 
manifold with dependence on /i and det {Taf^^ < 1 and let the non- 
degenerate homoclinic tangency be crossed at non-zero speed at fi = hq. 
Then for We > 0, 3(;Lti, C {/iq, iiq + e) and a residual subset J C (/ii, fi2) 
such that for fi ^ J , has infinitely many sinks. 

Models of such diffeomorphisms were constructed by Gambaudo and 
Tresser[ p^] and Wang proved that the Newhouse set has positive Hausdorff 
measure 
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After these results, intense research followed on the unfolding of homo- 
clinic tangencies and an essential question was whether, in addition to in- 
finitely many sinks, there would also be infinitely many strange attractors 
near the homoclinic tangencies. The question was positively answered by 



Colli [Q. The main result is: 

Theorem (Colli) Let /o G Dif f°°{M) be such that /o has a homoclinic 
tangency between the stable and unstable manifolds of a dissipative hyper- 
bolic saddle po. Then, there is an open set Q C Dif f°°{M) such that 

(a) /o G n 

(b) there is a dense subset D G Q such that for all / G -D , / exhibits 
infinitely many coexisting Henon-like strange attractors. 

Having established the existence of infinitely many sinks and infinitely 
many strange attractors near homoclinic tangencies, a question of practi- 
cal importance is the stability of the phenomenon under small random per- 
turbations of the deterministic dynamics. It turns out that the answer to 
this question is negative. Therefore under small random perturbations only 
finitely many physical measures will remain. 

Theorem (Araujo|2^) Let f : M — > M be a diffeomorphism of class 
C' , r > 1 , of a compact connected boundaryless manifold M of finite dimen- 
sion. If / = /a is a member of a parametric family under parametric noise 
of level £ > 0, that satisfies the hypothesis: 

There are K & N and i^o > such that, for all k > K and a; G M 

(A) fix.A)DB''ix),^o); 

(B) f^{x,p'^) « m ; 

then there is a finite number of probability measures ^i, ■■■ fii in M with 
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the properties 

1. ^i, - ■ ■ fii are physical absolutely continuous probability measures; 

2. suppz/jflsuppiUj for all 1 < « < J < / ; 

3. for all X G M there are open sets Vi = Vi{x) , ■ ■ ■ ,Vi = Vi{x) C A such 
that 

(a) VinVj = (lS, 1 < i < j <l ; 

(b) {A\{ViU---UVi)) = 0; 

(c) for all 1 < i < / and iy°°— a.e. t G we have 

1 n-l . 

for every (f) E C (M, R). Moreover the sets Vi(x), ■ ■ ■ , VJ(x) depend continu- 
ously on a: G M with respect to the distance d^, {A, B) = (AAB) between 
jyoo _ niodO subsets of A. 



3 Small dissipative perturbations of conser- 
vative systems 

Conservative systems have a large number of coexisting invariant sets, namely 
periodic orbits, invariant tori and cantori. By adding a small amount of 
dissipation to a conservative system one finds that some of the invariant sets 
become attractors. Of course, not all invariant sets of the conservative system 
will survive when the dissipation is added. However, for sufficiently small 
dissipation many attractors (mainly periodic orbits) have been observed in 



typical systems. Poon, Grebogi, Feudel, Hunt and Yorke p3 [Hi have 
extensively studied these effects in the single and double rotor, the Henon 
map and the optical cavity map. They find a large number of attractors for 
a small amount of dissipation, in particular in the double rotor map. The 
large number of coexisting stable periodic orbits has a complex interwoven 
basin of attraction structure, with the basin boundaries permeating most of 
the state space. The chaotic component of the dynamics is in the chaotic 
saddles embedded in the basin boundary. The systems are also found to be 
highly sensitive to small amounts of noise. The authors have argued that 
the two attributes, namely, accessibility to many different states and high 
sensitivity, are an asset in the sense that they are suitable for an easy control 
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of the complex system. The transition between different stable states, poses 
however delicate problems in view of the large chaotic transients in the basin 
boundary. The problem of migration between attractors and their stability 
in multiple-attractor systems has also been studied by other authorsp9[ [pO|] . 

All this work is very interesting, however most of results are based on 
numerical evidence. It would be desirable to have some control of the effects 
by rigorous mathematical methods. The techniques of deformation stability, 
to be discussed in the next section, might provide such a tool, in some cases 
at least. 



4 Deformation stability 

The basic idea is that, when dissipation is added to a conservative system, it 
is only a part of the total phase space that is related to the invariant sets of 
the dissipative system. The question therefore is to find the subsets in the 
conservative phase space that correspond to the stable sets in the dissipative 
system. In a sense the situation is similar to KAM theory and to the role 
played by the constants of motion in this theory. In KAM theory, the first 
integrals Iq defined everywhere in phase space for the integrable system, are 
deformed into a set of constants of motion 1^ which are defined only over 
a subset of sufficiently irrational tori. By analogy, a similar deformation 
stability may exist for dissipative perturbations of conservative systems, the 
stable domain of the constants of motion being, for example, the closure of 
a family of periodic orbits. In the papers where these ideas 



were developed, we were concerned with permanence under deformation of 
invariant sets, not necessarily attracting invariant sets. Most results may 
however be adapted to the search for attractors. As an example a result will 
be proved which gives a rigorous criterium for the existence of stable periodic 
orbits in a class of maps. 

Let an e— family of maps be 

x' = bx + y + f{x,y,e) 
y = y + 9ix,y,6) 

with /(x, y, 0) = g{x, y,0) = For e = the map has marginally stable 
periodic orbits of all periods. Under perturbation some of the orbits become 
stable ones. 
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Theorem If / and g are jointly in {x, y, e) with /(x, y, 0) = g{x, y, 0) = 
0, there is an e such that for \e\ < \e\ an interior orbit of period p of the un- 
perturbed map becomes a stable orbit of the perturbed map if and only 
if: 

il)Edegix(^\y^'\0) U= 

n=0 

(2) 

edej: {dMx^:\y^'\e) + (1 - b)dyg{x^^\y^'\e)} |,=o< 

n=0 

Proof: 

Iterating the map p times the orbit condition is: 

{lf-l)x + y'-^ + Up{x,y,e)-Pp{b) = 
Vp{x,y,e) =0 

with 

Upi^^y^^) = f^'^~^f{Xn,yn,e) 
n=0 

+ (1 - b)-' Eil - bP-''-')g{xn, yn, e) (2) 

n=0 

p-1 

Vp{x,y,e) = E g{xn,yn,£) 

n=0 

and (3p{b) is a polynomial in b. 

b < 1 and the condition imply, by the implicit function theorem, the 
existence of a solution x = h{y,e) of the first equation in (|1|). The second 
equation becomes 

Vp{h{y,e),y,e) = ev^{x,y,e) =0 (3) 

being a function. Condition (1) in the theorem follows from Eq.(^ for 
e = and the non- vanishing of the expression in the condition is required by 
the application of the implicit function theorem to (0). 

The eigenvalues of the p-iterated unperturbed map are 1 and bP. The sign 
in condition (2) is required to perturb the first eigenvalue to a smaller value. 
□ 
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5 Coupled quadratic maps 



A system as simple as one composed of two coupled quadratic maps, may 
have an infinite number of stable periodic orbits. Let the system be 

Xi{t + 1) = 1 - ((1 - c)xi(t) + CX2{t)f , s 

X2it + 1) = 1 - /i, (cXi(t) + (1 - C)x2it) f ^ ^ 

with X G [—1, 1], and = 1.401155... , which is the parameter value of the 
period doubling accumulation point. 

Theorem |^ For sufficiently small c there is an such that the system 



@ has stable periodic orbits of all periods 2" for n > N. 
Proof: 

Two essential features in the proof are the permanence of the unstabilized 
orbits in a flip bifurcation and the contraction effect introduced by the convex 
coupling. Only a sketch of the proof will be presented. For more details refer 
to Ref.[||. 

The bifurcations leading to the Feigenbaum accumulation point at /i* 
are flip bifurcations. This means that, after each bifurcation, the orbit that 
looses stability remains as an unstable periodic orbit. Therefore, (for c = 0) 
at /i = /i* the system has an infinite number of unstable periodic orbits 
of all periods p = 2". 

The proof has two basic steps. First one proves that, for sufficiently 
small c 7^ 0, these periodic orbits still exist in the system (|^). Second, that 
for any such c, there is an N such that there is at least one stable orbit 
for all periods jo = 2" with n > N. For both steps an important role is 

played by the instability factor, given by (^Z'-^-') (xp) at the fixed points Xp 
of the iterated map. Using the properties of the Feigenbaum - Cvitanovic 

functional equation one finds that the instability factor (^Z*^^-*^ (x*) converges 
to a fixed non-zero uniformly bounded value for all orbits. One now proceeds 
to the proof of the theorem. 

First step: Permanence of the periodic orbits for small c 

Let Xp G [—1, 1] X [—1, 1] be, for example, the coordinate of the p— periodic 

orbit closest to zero. The sequence {x*} = |x* : p = 2, 4, 8, ...| is an element 
of a £oo Banach space (sup norm). The collection of fixed point equations 

r{p) f * \ * ri 

Jfl* {^p: Cj Xp U 
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defines a C°°— mapping F{x*, c) from R ^ (.^ . Because (/^^^^ (x*) at 
c = is negative and bounded for all p, the derivative DiF of the mapping in 
the first argument is invertible. Therefore, by the implicit function theorem 
for Banach spaces, there is a c* such that for c < c* the function x*{c) : R — *■ 
(.00 is defined, that is, there are p— periodic orbits for all periods p = 2". 

For the uncoupled case (c = 0) the instability factor {^f'^^^^ {x*) for each 
mapping is the product (— 2//*)^ nfc=i over the orbit coordinates. For 
c < c*, the orbit structure being preserved, their projections on the axis are 
continuous deformations of the c = case which will preserve the geometric 
relations of the Feigenbaum accumulation point. Hence the same products 
for the projected coordinates suffer changes of order a{c)pX^ and remain 
bounded. 

Second step: Stabilization of at least one orbit for all periods p = 2" with 
n > N{c) 

The stability of the periodic orbits is controlled by the eigenvalues of the 

Jacobian Jp = ^^j^*'" in the fixed point of f^]c- The map (|) is a composition 
of two maps /i o 

■ 1 - 
1 - 





/i 

r f \ _^ ( {l-c)Xi+CX2 
^ \ ^2 J \^ CXi + (1 - c)Xi 

and by the chain rule the Jacobian has determinant 

det Jp = (1 - 2c)P(-2/x,)2f f[ x^{k)x2{k) (5) 



k=l 



Because of the permanence of the periodic orbits, for small c, in the neighbor- 
hood of the original coordinates (those for c = 0), the product of the last two 
factors in is uniformly bounded for all p. Then for all sufficiently large 
p, I det Jp\ < 1. The question is how this overall contraction is distributed 
among the two eigenvalues of Jp. 

To discuss the nature of the eigenvalues we may use a first order approx- 
imation in c. For a periodic orbit of period p = 2" we define 

(.) _ r (-2/i,)'^-'+i ui=i xm if g > n . , 

1 if/>o 
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For small c consider the linear approximation to the Jacobian Jp 

The eigenvalues are 
A± = i(l-pc)(xg + xg) 

(7) 

If the periodic orbit runs with the two coordinates xi and X2 synchronized 
then 

^+ = ^i,p 

A_ = (l-2pc)Xi,p 

and the orbit being unstable for c = it remains unstable for c 7^ 0. However 
if the two coordinates are out of phase by | steps the radical in Eq.(|^ is 



x{i)x{i + 1) ■ ■ ■ x{i H 1) 

i=i 2 

The existence of a superstable orbit for all periods jo = 2" implies that 
a.t n = fi^ the product 11^=1 -f^ik) has an odd number of negative-valued 
coordinates. Therefore the two eigenvalues are complex conjugate and, for 
small c, the contraction implicit in (j^) is equally distributed by the two 
eigenvalues. Therefore for sufficiently large N all orbits of this type with 
p > 2^ become stable periodic orbits. □ 

The attracting periodic orbits of the coupled system being associated to 
the unstable periodic orbits of the Feigenbaum cascade, the basins of attrac- 
tion will be controlled by the neighborhoods of these orbits, in each coordi- 
nate. Therefore a checkerboard-type structure is expected for the basins of 
attraction. 

References 

[1] S. Schiff et al.; Nature 370, 615 (1994) 



n 



k=l 



9 



[2] J. Foss et al.; Phys. Rev. Lett. 76, 708 (1996) 

[3] H. S. Seung; Neural Networks 11, 1253 (1998) 

[4] W. J. Freeman; Biol. Cyb. 56, 139 (1987) 

[5] J. A. Dente and R. Vilela Mendes; Phys. Lett A211, 87 (1996) 

[6] R. S. MacKay and J. - A. Sepulchre; Multistahility in networks of weakly 
coupled bistable units, Warwick preprint 44/1994. 

[7] H. Bohr, P. McGuire, C. Pershing and J. Rafelski; Threshold noise as a 
source of volatility in random synchronous asymmetric neural networks, 
Univ. Arizona report, Tucson 1997. 

[8] S. M. Hammel et al; J. Opt. Soc. Am. B2, 552 (1985) 

[9] M. Brambilla et al.; Phys. Rev. A43, 5114 (1991) 
[10] P. Marmillot et al.; J. Chem . Phys. 95, 1206 (1991) 
[11] J. P. Laplante and T. Erneux; Physica A188, 89 (1992) 
[12] K. L. C. Hunt et al.; J. Chem. Phys. 96, 7019 (1992) 
[13] F. Prengel et al; Phys. Rev. B50, 1705 (1994) 

[14] S. Yoden; Nonhnear Analysis, Methods and Apphcations, 30, 4607 
(1997) 

[15] E. Ott, C. Grebogi and J. A. Yorke; Phys. Rev. Lett. 64, 1196 (1990) 

[16] Focus Issue on Control and Synchronization of Chaos] Chaos 7, no. 4 
(1997) 

[17] S. Newhouse; Proc. A. M. S. Symp. Pure Mat. 14, 191 (1970) 

[18] S. Newhouse; Topology 13, 9 (1974) 

[19] S. Newhouse; Publ. Math. I.H.E.S. 50, 101 (1979) 

[20] J. Pahs and M. Viana; Annals of Math. 140, 207 (1994) 



10 



[21] C. Robinson; Comm. Math Phys. 90, 433 (1983) 

[22] J. M. Gambaudo and C. Tresser; J. Stat. Phys. 32, 455 (1983) 

[23] X.-J. Wang; Comm. Math. Phys. 131, 317 (1990) 

[24] E. ColU; Coexistence of inifinitely many strange attractors, IMPA report 
F-086/96, to appear in Ann. de 1' Institut Henri Poincare (Analyse non- 
hneaire) . 

[25] V. Arailjo; Attractors and time averages for random maps, IMPA report 
F-109/98. 

[26] L. Poon and C. Grebogi; Phys. Rev. Lett. 75, 4023 (1995) 

[27] U. Feudel, C. Grebogi, B. R. Hunt and J. A. Yorke; Phys. Rev. E54, 71 
(1996) 

[28] U. Feudel and C. Grebogi; Chaos 7, 597 (1997) 

[29] R. Weigel and E. Atlee Jackson; Int. J. of Bifurcation and Chaos 8, 173 
(1998) 

[30] K. Kaneko; Phys. Rev. Lett. 78, 2736 (1997) 

[31] R. Vilela Mendes and J. T. Duarte; J. Math. Phys.; 24, 1772 (1983) 

[32] R. Vilela Mendes and J. T. Duarte; Lett. Math. Phys. 6, 249 (1982) 

[33] R. Lima and R. Vilela Mendes; Nonlinear Analysis, 12, 1061 (1988) 

[34] R. Vilela Mendes; in Deformation Theory of Algebras and Structures 
and applications , pp. 981 - 1013, M. Hazewinkel and M. Gerstenhaber 
(Eds.) Kluwer 1988 

[35] R. Carvalho, R. Vilela Mendes and J. Seixas; Physica D126, 27 (1999) 



11 



